Granular impact and the critical packing state 
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Impact dynamics during collisions of spheres with granular media reveal a pronounced and non- 
trivial dependence on volume fraction 0. Post impact crater morphology identifies the critical packing 
state (j)cps^ where sheared grains neither dilate nor consolidate, and indicates an associated change 
in spatial response. Current phenomenological models fail to capture the observed impact force 
for most 0; only near (j)cps is force separable into additive terms linear in depth and quadratic in 
velocity. At fixed depth the quadratic drag coefficient decreases (increases) with depth for cj) < (j)cps 
(0 > (t)cps)- At fixed low velocity, depth dependence of force shows a Janssen-type exponential 
response with a length scale that decreases with increasing cj) and is nearly constant for > (t)cps- 

PACS numbers: 96.15.Qr,96.20,45.70.-n,83.80.Fg,47.50.-d 
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Impact of objects into unconsolidated granular mate- 
rials [l| like sand is relevant in many settings (e.g. ter- 
minal ballistics of projectiles, hammered intruders, and 
ground- foot interaction \3^) and is of scientific interest be- 
cause the localized strain field around the impactor gen- 
erates interaction between fluid and solid granular states. 
Unlike rapid granular flows which can be described by 
hydrodynamic-like equations [3] (e.g. steady chute flow 
with no enduring contact networks), equivalent compre- 
hensive and tested continuum descriptions for the mixed 
solid/fluid regime are lacking. 

In the absence of governing equations and with a 
scarcity of direct force measurements, numerous phe- 
nomenological force models have been proposed over hun- 
dreds of years [J-[l4 ^^ explain the observed dependence 
of penetration depth, crater morphology, and collision 
duration on impact velocity and intruder and grain prop- 
erties. To the best of our knowledge, all existing models 
assume that the granular resistance force F can be sepa- 
rated into independent functions of position and velocity 
such that F{z^v) = Fz{z) + Fy{v)^ where z and v are 
the projectile's depth below the initial free surface and 
velocity respectively. The depth dependent term Fz has 
been modeled as constant Jj,l8|, [ll], linear [5, 8, 14] and 
as a modified exponential [l2|. Velocity dependence also 
remains uncertain; it has typically been treated as an in- 
ertial drag, Fy = av'^ i,i,|lll,[l3,[l3,[l5l, although linear 
^8] and constant [16l] forms have also been proposed. 

No impact experiments or models have systematically 
examined the effect of volume fraction (/), a parameter 
that largely determines the mechanical response of slowly 
sheared granular media [17]. (/>, the ratio of material vol- 
ume (total mass divided by constituent solid density) to 
occupied volume, ranges between 0.55 and 0.64 [18] for 
dry non-cohesive granular media of slightly polydisperse 
and nominally spherical grains. In slow granular shear- 
ing, ^ is a principle determinant of yield stress through 




FIG. 1: (Color online) Influence of volume fraction on crater- 
ing at vo = 257 ± 3 cm/s. (a) Impact crater at = 0.61 and 
laser line, (b) Surface displacement h relative to grain diame- 
ter d increases with volume fraction (0 = 0.579, 0.589, 0.600, 
0.610, and 0.622). (c) Post-impact change in bed volume Ay 



relative to sphere volume Vs vs. is at 



= 0.591, indi- 



cating the location of the critical packing state. 



changes in flow structure: for large ^, materials dilate and 
flow locally in shear bands while for small cj) they consol- 
idate and flow globally without shear bands [17]. Only 
at the critical packing state [19] or CPS, which occurs 
at intermediate (/), is the volume fraction constant under 
shear. Volume fraction effects are expected for penetra- 
tion as well, and recent work examining slow constant 
velocity intrusion of a cylinder has observed signatures 
of a phase transition at (/) ~ 0.60 [20|. 

Here we directly measure the time resolved impact 



force and post-impact crater profile to determine the in- 
fluence of <j). In our experiments an R = 1.98 cm ra- 
dius steel sphere with total mass m = 270 g is dropped 
with initial collision velocities < vq < 350 cm/s into 
the center of a 24 x 24 cm^ cross-section box filled to a 
depth of ~ 30 cm with d = 300 jam mean diameter glass 
spheres (density p = 2.52 g/cm^). An accelerometer on 



the sphere records vertical acceleration a (see Ref. [14 1 
for details) from which the non-dimensional penetration 
force is determined: F = F/mg = a/g + 1, where g is 
the acceleration due to gravity. The sphere's velocity is 
v{t) = '^0 + /o ciit')dt' and its lowest point beneath the 
initial sand surface is z{t) = Jq v(t')dt' ^ where t = 
corresponds to the time the bottom of the sphere first 
contacts the bed. Impactor and glass sphere dimensions 
ensure that our results are not influenced by finite size 
effects 17| or interstitial air [21]. Bed dimensions were 
chosen to eliminate boundary effects [22|, |23| , but as we 
will suggest, the influence of horizontal walls appears to 
increase with increasing (j). 

States with 0.57 < (j) < 0.63 are generated by initially 
flowing air upward through the rigid distributor base of 
the bed to create a fully fluidized state, then decreasing 
air flow below fluidization onset and vibrating the con- 
tainer to reduce (j) to the desired value J2|. Air flow and 
vibration are stopped during bed height measurement 
and after reaching the desired volume fraction: collisions 
occur in a quiescent bed. (j) is measured to a precision 
of 0.001 using an ultrasonic range finder to determine 
bed height. Additional x-ray absorption measurements 
confirmed that vertical variation in (j) is small (less than 
0.004 for A^ = 0.04) as in previous observations, e.g. 
|24| . Surface profiles are measured using laser line pro- 
filometry [25|, see Fig. [TJa). 

We first examine the impact crater (studied previously 
only at fixed 0, e.g. [26]) for constant vq and show that 
its shape and the amount of displaced material are sen- 
sitive to volume fraction. Figure [IJb) indicates that at 
large (j) the crater has a high rim and a small central peak 
(the remnant of a granular jet [9]), while at small (j) the 
crater is deeper and has a lower rim and a larger cen- 
tral peak. To quantify these changes, we use the crater 
height profile /i(r), where r is the radial distance from 
the crater center, to calculate the post-impact change in 
bed volume /SV = 27r Jq rh{r)dr — Vs vs. 0, where Vg is 
the sphere volume (the sphere is always fully submerged 
post-impact). Figure [TJc) shows that /SV is negative at 
low (j) and positive with smaller slope at high (j). /SV can 
be interpreted as the volume of grains disturbed in the 
collision, V^, times the average change in volume fraction 
within Vd. Since V^ > 0, the transition in /SV from nega- 
tive (compaction) to positive (dilation) indicates that the 
critical packing state (no change in (j) with shear) occurs 
at a volume fraction of (j)cps = 0.591 ± 0.005 for the glass 
beads used in this study. We quantify the proximity to 
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FIG. 2: (Color online) Effects of volume fraction on impact. 
(a) Penetration depth decreases with increasing Acj) while (b) 
collision duration is nearly unchanged [vq = 134 cm/s (blue 
circle) and vo = 283 cm/s (red triangle)]. In (a) and (b) 
changes in response occur near A(/) = 0. (c) Penetration force 
increases with time for A(/) = —0.007 (green solid) but de- 
creases for Ac/) = +0.015 (orange dashed) at vq = 190 cm/s. 
Curves are averages of 10 experiments; gray regions are zblcr 
and indicate the high degree of experimental repeatability. 
Inset: zv-p\ane trajectories differ with A0 for the same vq. 



GPS with A(j) = (j)- (t)cps' 

At the critical packing state the granular medium is 
effectively incompressible, suggesting the possibility of 
simpler impact dynamics in its vicinity and qualitative 
differences in dynamics between states with (j) < (j)cps 



and 



> 



due to (j) dependent changes in flow as in 



Fig. life). Beginning with kinematics. Fig. [2fa) shows 
that the penetration depth at fixed vq decreases with in- 
creasing volume fraction as expected. Penetration depth 
decreases sharply for Acj) < but more gradually for 
A(j) > with a total decrease of ~ 30%. Collision du- 
ration tc, although it varies little, is more sensitive to 
variation in cf) for A^ < 0, see Fig. EJb). A^ influences 
kinematics over the course of the collision as indicated by 
impactor trajectories in the zv-pldiiie [inset of Fig. [21(c)]. 

The (f) dependence of kinematics originates in the 
impact force, which, as Fig. [2](c) shows, differs above 
and below the critical packing. For example, at vq = 
190 cm/s and below GPS (A^ = -0.007), F increases 
with time from 2.5 to 4 times the weight of the in- 
truder (like slow penetration in granular media |l7J), 
while above GPS (Acj) = +0.015), F decreases by about 
the same amount over the same time interval. In both, 
collision onset is marked by an initial jump in F, and col- 
lision termination is characterized by a sharp decrease in 
force of magnitude Fgtop [14]. The average slope of F{t) 
during collision decreases with increasing (j) (and vq^ see 
[ij]); the decrease is more rapid for A(/) < 0. 

To better characterize the dependence of impact force 
on volume fraction, depth, and velocity, we measured 
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FIG. 3: (Color online) Isoforce contours (boundaries of col- 
ored regions) and impactor trajectories (circles, 6 ms sam- 
pling interval, VQ/\fRg = 2, 4, 6) in the z'^-plane differ 
(a) below, (b) at, and (c) above the critical packing state. 
At CPS, contours predicted by Eqn. [T] (white curves) best 
match experiment, and /\(j) dependent deviations from Eqn.[Tl 
AF = F — Fu, are least localized. Eqn. [T] fit parameters: 
k = 0.62, 0.88, 1.21 and a = 0.76, 0.82, 1.07 for A0 = -0.01, 
0.00, 0.02, respectively. 



F{t) for —0.01 < A(j) < 0.03 varied in increments of 
0.002 and for < vq/ y/Rg < 6 varied in increments of 
« 5 cm/s. For each sampled time in each collision we 
calculated the position z{t) and velocity v{t) of the im- 
pactor; note that zi'-trajectories for distinct Acj) and vq 
do not intersect, see Figs.[2fc) and[3l We then partitioned 
the Z'U-plane into 0.2 cm x 5 cm/s regions and calculated 
the average force in each region to find F{z,v, ^0). Fig- 
ure [3] presents isoforce contours of F{z^v, A^) for A^ 
below, at, and above CPS and reveals A^ dependent 
changes in F. For small z, force increases more rapidly 
with V as A0 is increased, while for small v, F increases 
more slowly with depth for increasing A^. Impactor tra- 
jectories (circles in Fig. [3]) show that with increasing 
depth F generally increases for Acj) < 0, decreases for 
A(/) > 0, and changes least at A^ = 0. 
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FIG. 4: (Color online) Variation of a^ with depth, (a) F 
and fits of F oc v^ (curves, see text) for Acj) below, at, and 
above CPS at fixed z/R = 1.1. (b) Variation of a with depth 
depends on A0. Error bars increase as the range of measured v 
decreases with z. Only near Acj) = is a^ depth independent. 
Error bars in (b) are 95% confidence intervals. 



We compare our data to a model from a recent study 
by Katsuragi et al. [15] of sphere impact into glass beads 
with initial (j) = 0.590. Using force data derived from 
high-speed imaging of impactor position vs. time, the 
authors proposed a "unified" force law of the form 






(1) 



where k and a are dimensionless constants. We fit k 
and a to F{z,v) at each A^ and found that Eqn. 1 
best describes F near CPS, see Fig. [3l Below and above 
CPS, differences between isoforce contours of experimen- 
tal data (boundaries between colored regions) and the 

model (white curves), v/^/Rg = y {Fu — ■^z)/a, are 
greatest at low v. Insets in Fig.[3]plot differences between 
the experiment and model fits, AF = F — Fu-, and show 
that for A^ < 0, the model underestimates F at both 
small z and small v, while for A^ > 0, it overestimates 
F at large z and small v. 

Figure [3] suggests that a force model linear in depth 
and quadratic in velocity is insufficient away from CPS. 
Isolating the velocity dependent contribution by examin- 
ing data at fixed depths (e.g. along vertical lines in Fig.[3j) 



using F = C 



Rg 



? , where ol and C are constants free 



to vary with depth, shows this is the case (see Fig. HJ^a) 
for example fits at zjR = 1.1). Figure IH^b) shows that 
(J decreases with z for A^ < and increases with z for 
A(\) > 0, which rules out a separable penetration force 
with a purely v'^ velocity dependence. Only near CPS is 
a^ independent of depth as in most granular impact-force 
models J4l4l4j. We note, however, that by adding a term 
linear in velocity, i.e. F = C -\- —^v + W^'^^ ^ becomes 
always positive and nearly independent of depth at all 0. 
The effects of proximity to CPS are also evident at low 
velocity where depth dependent frictional forces are ex- 
pected to dominate. F^^op (see Fig. [2]) characterizes the 
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FIG. 5: (Color online) /\(j) driven changes in grain force at the 
end of collision, (a) Fstop is linear in depth at the lowest /\(j) 
but becomes increasingly sub-linear as A(/) increases; curves 
are fits of Fstop vs. z/ R to a two parameter model similar to 
Janssen's Law (see text), (b) Characteristic length A from 
model decreases rapidly with /\(j) below CPS but is nearly 
constant above it (dashed curve is an exponential fit). Error 
bars in (b) are 95% confidence intervals. 



low velocity response at the end of collision and is plot- 
ted in Fig. [5ja) vs. depth for representative A(/). For the 
smallest A^, Fstop increases linearly with depth while at 
larger A^ the dependence is sub-linear. The force in the 
depth dominated regime for all A0 is well modeled by 
a single Janssen-like function J27| Fgtop = k' (l — e~^/^) 
[black curves in Fig. [5fa)], where A is a characteristic 
length and k' a constant. A decreases rapidly from ^ 30i? 
to ^ 2R as A^ approaches zero from below, see Fig.[5fb). 
For A^ > 0, A is nearly constant. In the limit of large A 
(low A0), the model gives a linear response Fstop ^ x 
with k = k' /\ as in Eqn.[Tl while in the opposite extreme, 
Fstop is constant as in Ref. [Uj . To determine if A is also 
dependent on container size, we reduced the width of the 
fluidized bed from 24 cm to 12 cm and found that A de- 
creased by ^ 2/3 at all A^. Our findings suggest that 
as (j) is increased, grains exert increasingly larger forces 
on the side walls, and presumably also on the intruder, 
at shallower depths which ultimately reduce the gravita- 
tional forces on the grains leading to a net reduction in 
the depth dependent component of the penetration force. 
Whether or not Fstop is ultimately linear in z at large A0 
in an unbounded container is an open question. 

By varying the volume fraction, we have shown that 
the dynamics of impact in granular media are richer than 
previously thought and that existing separable models of 
granular impact linear in depth and quadratic in veloc- 
ity do not capture all the details, likely due to changes 
in flow and the influence of boundaries associated with 
compaction and dilation. Dynamics in spatially extended 
systems based on ODE's are often incomplete; in the case 
of granular impact, they fail to explain why and how 



heuristic parameters such as k and a change with depth, 
velocity, and (j). To completely characterize granular im- 
pact it is likely that the full spatio-temporal response of 
the granular ensemble is needed; continuum approaches 
capturing the transition from static to flowing regimes 
appear promising, see for example [28|. 
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